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PARAMETERIZATION  OF  THE  ATMOSPHERIC  HEATING  RATE 
FROM  15  TO  120  KM  DUE  TO  Oa  AND  03 
ABSORPTION  OF  SOLAR  RADIATION 

INTRODUCTION 

In  order  to  numerically  simulate  the  earth's  upper  atmospheric 

circulation,  accurate  and  computationally  efficient  parameterizations 

of  solar  radiative  heating  are  required.  Of  major  interest  is  the  solar 

radiative  heating  by  ozone  absorption  in  the  Hartley  region  (2000- 

3000$) , Huggins  bands  (3000-3500&)  and  the  Chappius  bands  (4500-7500A) . 

Above  the  mesopause  0^  absorption  of  solar  uv  radiation  in  the 

Schumann- Runge  bands  (1750-2050A)  and  Schumann- Runge  continuum 
o 

(1250-1750A)  must  also  be  included. 

For  ozone  heating  rates  two  different  parameterizations  have 
been  previously  developed.  Lindzen  and  Will  (1973)  assumed  that  the 
0^  absorption  cross  sections  are  either  constant  or  have  an  exponential 
variation  over  individual  wavelength  intervals.  With  these  assumptions, 
they  obtained  simple  analytic  formulas  for  the  0^  heating  rate.  The 
advantage  of  this  method  is  that  improved  solar  flux  and/or  cross 
section  data  can  easily  be  Incorporated  into  the  parameterization 
scheme.  An  alternate  approach  has  been  developed  by  Lads  and  Hansen 
(1974)  for  general  circulation  models.  They  found  analytic  expressions 
for  the  0^  heating  rate  profiles  by  numerical  fits  to  exact  profiles. 
However,  to  incorporate  revised  solar  flux  and/or  cross  section  data 
this  method  requires  a repetition  of  their  fitting  procedures  to 
obtain  new  numeric*.!  coefficients.  This  approach  is  obviously  not  as 
flexible  in  some  res.  . cs  as  the  Lindzen  and  Will  (1973)  formulas. 
Accordingly,  the  lact<*  is  adopted  to  represent  0 2 and  0^  heating  rates 
by  Insolation  from  15  to  120  km. 


Note:  Manuscript  submitted  October  18,  1976. 
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EXACT  HEATING  RATE  CALCULATION 


To  develop  parameterlzatlons  for  any  atmospheric  circulation 
model,  it  is  best  to  have  "exact"  or  detailed  calculations  available  to 
assess  their  accuracy.  The  solar  fluxes  and  the  02  and  0^  absorption 
cross  sections  adopted  for  the  exact  calculation  are  tabulated  in 
Table  1.  Since  we  will  indicate  how  to  incorporate  improved  fluxes  and 
cross  sections  in  the  parameterlzatlons,  we  do  not  need  to  justify  the 
selected  input  data  (e.g.,  the  controversial  solar  flux  values  in  the 
1750-2050X  region). 

The  solar  fluxes  were  selected  from  the  following  sources: 
Broadfoot  (1972),  Brueckner  et  al.  (1976),  CIAP  Monograph  1,  Donnelly 
and  Pope  (1973)  and  Smith  and  Gottlieb  (1975).  The  02  and  0^  cross 
sections  are  based  on  Blake  (1973),  CIAP  Monograph  1,  Hudson  (1971)  and 
Hudson  and  Mahle  (1972).  Calculation  of  the  02  and  03  heating  rates  is 
straightforward  for  all  wavelengths  of  interest  with  the  exception  of 
solar  radiation  penetration  through  the  absorbing  0^  Schumann- Runge 
bands  (1750-2050A) . McConnell  (1974)  has  noted  that  transmission  of 
solar  radiation  in  each  band  i_,  Tr^  can  be  adequately  represented  by 


Tr1  - exp  |”(YiN2  + 


where  N2  is  the  column  density  of  0^  along  the  solar  radiation  path, 
and  5^  are  coefficients  given  in  Table  2.  The  total  0,  heating  rate 
in  the  Schumann-Runge  bands  (SRB)  is 


QSRB(02)  " £ Fi  0i  Tri  [°2J 


where  is  the  solar  flux  integrated  over  each  band,  [(^l  is  the  C>2 
number  density  and  the  cross  section  (o^)  is  given  by 


“i  + 


with  and  6.^  tabulated  in  Table  2.  An  accuracy  of  within  20  percent 
of  the  experimental  data  of  Hudson  and  Mahle  (1972)  can  be  achieved  with 
these  expressions  and  coefficients.  The  O2  column  density  was  computed 
using  a 16  point  Gaussian  quadrature  integration  scheme. 

Net  heating  rates  (eQ)  defined  as  the  energy  absorbed  minus  the 
chemical  energy  stored  by  the  dissociated  products  were  computed  with 
efficiency  factors. 


„ / he  _ \ ,hc 

V's 


where  D is  the  dissociation  energy  of  the  absorbing  molecules  (1.05  ev 
for  0^  and  5.12  ev  for  O2) , Xj  is  the  wavelength  of  the  photon,  h is 
Planck's  constant  and  c is  the  speed  of  light. 
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PARAMETERIZATION  OF  THE  03  HEATING  RATE 


For  the  Chappius  bands  (c)  we  follow  Lindzen  and  Will  (1973) 
and  write  the  0^  heating  rate  (Qc)  as 


ic_ 

[o3] 


F a 
c c 


-acN3 


(5) 


where  [03J  is  the  03  number  density  and  N3  its  column  density  along  the 
solar  radiation  path.  We  find  that 


5 -2 

F * 3.7  x 10  ergs  cm  sec 
c 

a * 2.85  x 10“21  cm2 
c 


(6) 


gives  the  best  fit  to  exact  calculations  of  Q . O is  0.62  times  the 

o c c 

peak  03  cross  section  at  ~6000A  and  Fc  is  0.7  times  the  integrated 

flux  over  the  Chappius  bands.  The  net  heating  rate  e Q is  obtained 

9 c c 
with  Xj  =*  6000A  for  the  efficiency  factor  from  equation  (4) . 

A similar  approximation  is  appropriate  for  the  Hartley  region 
(Ha)  from  2425-2775A.  Thus 


FHa  aHa  8 


"aHAN3 


(7) 


where 


-2 

Fa  ■ 5460  ergs  cm  sec 
Ha 

a„  ■ 8.8  x 10  ^ cm2 
Ha 


(8) 
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with  a equal  to  0.8  times  the  peak  0„  cross  section  at  ~2525&  and 
Ha  J q 

FHa  the  integrated  flux  over  the  2425-2775A  interval.  The  net  heating 
rate  is  calculated  with  X^  » 2500&  in  equation  (4) . 

In  the  Huggins  bands  (Hu,  defined  here  as  the  2775-3600& 
interval)  we  assume,  as  Lindzen  and  Will  (1973),  that  the  0^  cross 
section  has  an  exponential  variation 


a **  a„  e 
hu 


-MX 


(9) 


To  obtain  greater  accuracy  than  Lindzen  and  Will  (1973),  we  break  the 
Huggins  bands  up  into  two  intervals.  The  heating  rate  is 


3600&  long 

f "CTN3  f 

/ Ix  a e J dX  + I z2  0 


xHu 

[o3J 


-CTN„ 


dX 


long 


short 


or 


^Hu  1 1 t . T ^longl 

toji ' 55;  J1! + 'yv  exp  l-w  J 


-I2 


-MX 


short 


]l 


(10) 


where  the  atmosphere  is  optically  thin  at  3600&.  Numerically,  an 
excellent  fit  to  the  exact  heating  rate  is  obtained  with 
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1 

2 
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X 


short 


59.2 

ergs 

-2  -1 
cm  sec 

tr1 

49.3 

ergs 

-2  -1 
cm  sec 

r1 

0.0127 

r1 

2805 

i 

0.0125  cm 


The  intensities  I..  and  I-  can  be  related  to  fluxes  integrated  over  the 

^ o 

respective  wavelength  intervals.  AA^  * 3600  - ^^ong  * 345A  and 


AA,  ■ A -A  * 25oX.  Thus,  F.  “ I AA,  ■ 2.04  x 10^  ergs 

2 long  short  * 1 1 1 ° 

“2  A -2  “1 

cm  sec  and  F-  * I-AA,  * 1.23  x 104  ergs  cm  sec  . From  Table  1 

1 L * 4-2-1 

the  integrated  flux  over  interval  1 is  4.4  x 10  ergs  cm  sec  and 

F^  is  0.46  times  this  value,  while  the  integrated  flux  over  the 

2775-3055A  interval  is  1.25  x 10^  ergs  cm  2sec  ^ and  F^  is  equal  to 

this  flux.  In  interval  1 the  atmosphere  is  partially  optically  thin 

and  all  photons  are  not  absorbed. 

In  the  Herzberg  continuum  (2060-2425^)  both  0^  and  0^  absorb 
solar  radiation  and  the  principal  region  of  heating  occurs  at  35-55  km. 
Adequate  representation  of  this  heating  (QHz)  can  be  obtained  with 

QHz  " FHz  {aHz(02)[02'  + aHz<03)^°3'}  eKp  ['0Hz(.02)K2 

-ffHz<°3»,3]  (l3 


°Hz<02) 


°az<03> 


3 -2  -1 
1.5  x 10  ergs  cm  sec 


in- 24  2 

6.6  x 10  cm 


4.9  x 10" 18  cm2 


■ 2290&  to  calculate 


where  F is  0.79  times  the  integrated  flux  in  this  wavelength  region. 
Hz 

The  parameterizations  for  the  Chappius,  Huggins,  and  Hartley 
bands  are  accurate  to  better  than  2 percent.  With  the  inclusion  of 
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the  Herzberg  continuum  parameterization  the  accuracy  for  the  summed 
C>2  heating  rates  decreases  to  5 percent,  primarily  because  equation  (12) 
underestimates  the  actual  heating  rate  below  42  km.  In  Figures  1 and 
2 the  individual  contributions  to  the  total  and  net  heating  rates, 
respectively,  are  illustrated  based  on  the  above  parameterizations . 

They  are  diurnally-averaged  heating  rates  at  the  equator  during  equinox. 
The  model  atmosphere  used  in  the  computations  is  given  in  Table  3. 
Conversion  factors  for  heating  rates  in  °K  day  ^ to  ergs  cm  ^sec  ^ are 
given  in  Table  4. 

Above  45  km  0^  absorption  in  the  Hartley  region  is  the  dominant 
heat  source,  whereas  below  28  km  the  Chappius  band  absorption  is 
dominant.  In  the  intermediate  region  the  Huggins  bands  dominate. 

The  secondary  peak  in  the  heating  rate  at  ~85  km  corresponds  to  the 
secondary  maximum  in  the  0^  concentration  in  the  model  atmosphere. 

Thus,  thermodynamically,  this  [0^]  peak  is  expected  to  be  quite 
important  and  its  absolute  value  must  be  accurately  known.  Absorption 
in  the  Herzberg  continuum  makes  a small  contribution  to  the  heating 


PARAMETERIZATION  OF  THE  02  HEATING  RATE 

The  Schumann-Runge  continuum  (SRC)  can  be  split  up  into  two 

o -17  2 

main  regions:  1250-1520A  where  the  CL  cross  section  —10  cm  , and 

Q ^ ""MX 

1520-1750A  where  the  02  cross  section  is  proportional  to  e . For  the 

1250-1520&  region  we  may  approximate  the  heating  as 


with 


XSRC 

[°2] 


F a e 
SRC  SRC 


aSRCN2 


„ , . -2-1 
fsrc  “ 1,1  erg  cm  sec 

aSRC  = 1 X 10’17  Cm2 
eSRC  = ^-41  ^or  net  Seating 


(13) 


where  F___  is  the  integrated  flux  over  this  wavelength  region  and  a 

SRC  bKJ 

is  the  average  cross  section.  For  the  1520-1750A  region  we  divide  it 

-MA 

into  two  intervals  with  cr(C>2)  « e • and  obtain  a total  heating  rate 
of 


XSRC 

[o2] 


1_ 

N„ 


~aZN2  . 1s~1SL  "amN2  Xs  *asN2 

_e  + M 6 're 


(14) 


where  denotes  long  wavelength  end  of  spectrum,  £ denotes  short  wave- 
length end  and  a denotes  CL  cross  section  at  the  junction  of  these 

m 2m 

intervals  (— 1660A).  The  following  values  were  found  to  give  an 
excellent  fit 
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a 

m 


_2 

3.43  ergs  cm  sec 


-2 

1.35  ergs  cm  sec 


-19  2 

2.9  x 10  cm 


i c/  in*18  2 

1.54  x 10  cm 


-1 


-1 


(15) 


a 

s 


1.1  x 10 


-17 


2 

cm 


— is  approximately  0.6  times  the  integrated  solar  flux  in  the 
™ _ £ 

1660-1750A  interval,  while  is  approximately  0.5  times  the  solar  flux 

in  the  1520-1660&  interval.  The  0„  cross  sections  a. , a , and  0 are 

2 l m s 

approximately  the  values  at  1750,  1660,  and  1520A,  respectively.  To 
compute  the  net  heating  rate  we  recommend 


el 

s 

M 


a 

m 


a 

s 


-2 

0.98  ergs  cm  sec 
-2 

0.43  ergs  cm  sec 

-19  2 

2.9  x 10  an 

. 7 in-18  2 
1.7  x 10  cm 

-17  2 

1.15  x 10  cm 


(16) 


where  £ cs  0.3. 

The  main  contribution  to  atmospheric  heating  from  the  0£ 
Schumann-Runge  bands  (SRB)  occurs  in  the  60-100  km  region.  In  view 
of  the  uncertainty  in  the  solar  flux  magnitude  in  the  Schumann-Runge 
bands  coupled  with  the  complexity  of  its  transmission  through  the 
atmosphere,  we  demand  only  an  accuracy  of  + 20  percent  in  this 
parameterization.  The  following  expression  represents  the  total 


1 


heating  rate  in  the  60-100  km  region 


SRB 


[02]  aN2  + bN^: 


where 


a **  0.143  cgs  units 

g 

b * 9.64  x 10  cgs  units 


and  if  N2  < 101®  cm^,  then 


SRB  Q ---9  -1 

9.03  x 10  ergs  sec 


[0,] 


For  the  net  heating  rate  we  obtain 


a = 0.67  cgs  units 


b = 3.44  x 10  cgs  units 


and  if  N < 10^  cm^,  then 


SRB  „ in-19  -1 

~ 2.43  x 10  ergs  sec 


[o2] 


(17) 


(18) 


(19) 


(20) 


(21) 


By  appropriate  scaling  of  both  coefficients,  heating  rates  for  alternate 
solar  flux  values  can  be  obtained. 


The  individual  contributions  to  the  02  heating  rates  from  the 
above  parameterizations  are  shown  in  Figures  1 and  2.  Near  60  km  the 
total  02  heating  rate  is  accurate  only  to  + 25  percent,  but  above  75  km 
it  improves  to  within  5 percent.  Heating  in  the  Schumann- Runge  bands 
is  the  dominant  heat  source  between  88  and  96  km,  while  absorption  in 
the  SR  continuum  is  most  important  above  96  km.  The  relative 
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1 


10 


MBPj pp 


importance  of  the  SR  bands  in  the  80-90  km  region  depends  significantly 
on  the  0^  concentration  there. 

The  overall  accuracy  of  the  parameterized  heating  rates  is 
+ 5 percent  over  the  altitude  region  15-120  km.  Below  approximately 

3 

80  km  the  0(  P)  formed  by  0^  and  0^  dissociation  recombines  quickly, 

and  the  actual  heating  rate  of  the  atmosphere  is  the  total  photon  energy 

3 

absorbed.  Above  80  km  the  0(  P)  produced  by  0^  and  0^  dissociation 
may  be  transported  significantly  in  the  vertical  direction  before 
recombining,  and  the  actual  heating  rate  is  the  net  heating  rate  (eQ) 
plus  the  chemical  energy  released  by  locally  recombining  0(  P) . To 
a good  first  approximation  we  can  estimate  the  actual  heating  rate 
by  summing  the  total  0^  heating  rate  and  the  net  0^,  heating  rate.  In 
Figure  3 this  estimated  actual  heating  rate  is  illustrated  for 
solstitial  conditions.  Only  the  latitudinal  variation  of  solar 
radiation  was  included;  the  model  atmosphere  was  invariant  with 
latitude.  Our  heating  rate  agrees  well  with  the  results  of  Park  and 
London  (1974)  when  the  comparison  is  made  at  the  same  density  levels. 
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PARAMETERIZATION  OF  DIFFUSE  SCATTERED  SOLAR  RADIATION 


Solar  radiation  longward  of  3000^  is  not  strongly  absorbed  in 
the  earth's  atmosphere  and  can  undergo  multiple  scattering  in  the 
atmosphere  by  molecules  and  particles.  Fortunately,  in  atmospheric 
heating  calculations  the  major  scattering  effects  are  in  the  Chappius 
bands  region  of  the  solar  spectrum  (Lacis  and  Hansen,  1974) . Their 
results  strongly  suggest  that  the  diffuse  solar  radiation  can  be  modeled 
by  a pure  0^  absorption  region  on  top  of  a reflecting  layer  with  an 
effective  albedo  (S^)  that  depends  on  the  ground  reflectivity  (R^)  and 
lower  atmosphere  albedo  (Rfl)  as  follows: 


0.856  R 

uS1  * R + [1  - R ] 

o a 1 aJ  1 - 0.144  R 


(22) 


g 


Let  the  vertical  optical  depth  due  to  0^  absorption  be 


O 

/ 


and  z 


(23) 


where  n is  the  0^  number  density,  a is  its  cross  section,  z is  height, 
and  dx  ■ -andz.  With  9 as  the  solar  zenith  angle  and  s as  the  distance 
along  the  solar  radiation  path  (increasing  in  positive  value  with 
decreasing  z),  we  define  y * cos9  and  yds  * -dz.  The  intensity  (I)  of 
the  solar  radiation  is  separated  into  upward  (I+)  and  downward  (I  ) 
components.  The  associated  flux  is 


F - 


(24) 
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for  a plane  parallel  atmosphere.  For  visible  sunlight  in  the  purely 
absorbing  atmosphere  the  equation  of  transfer  is  simply 


y dT  1 (26) 

since  atmospheric  thermal  emission  is  negligible.  Equation  (26)  has 
the  solution 

X 

I - I0  eU  . (27) 


For  the  downward  intensity  component  the  boundary  condition  is  applied 
at  the  top  of  the  atmosphere  and  is 

r<T  - 0,y)  - Fq  5(y  + yQ)  (28) 

where  F0  is  the  downward  directed  solar  flux  with  direction  y^. 

Thus 

2 

l'(T,y)  - F0  6(y  + yQ)  eU  , y < 0 . (29) 

For  the  upward  intensity  component  the  boundary  condition  is  applied 
at  the  reflecting  surface  (t  ■ T*)  and  is  a Lambert  surface  with 
albedo  (15* ) , 


represents  the  transmission  of  the  reflected  light  off  the  surface. 
This  second  term  in  equation  (33)  is  thus  the  heating  due  to  diffuse 
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reflected  and  scattered  solar  radiation,  whereas  the  first  term  in 
equation  (33)  represents  direct  solar  radiation  heating. 

Lacis  and  Hansen  (1974)  recommend  a simple  approximation  for 
the  transmission  function  E2.  It  is 


E2(t*-t)  =s  e 


-(t*-t)M 


where  M is  the  magnification  factor  for  the  vertical  optical  depth  and 
has  an  effective  value  of  1.9  for  diffuse  radiation. 

Comparisons  of  their  detailed  calculations  with  the  results 
obtained  from  the  following  expression 


T 

P 


+ 2<?0u0e 


0 -1.9(t*-t) 1 

6 1 


show  good  agreement  for  diffuse  radiative  heating  (within  10  percent) . 

If  the  effective  albedo  is  a3^  * 0.25,  a reasonable  global  average,  then 
the  0^  heating  rate  from  Chappius  bands  absorption  is  increased 
by  30  percent  throughout  the  ozone  layer  due  to  diffuse  solar  radiation. 
This  description  of  diffuse  radiation  should  also  prove  adequate  for 
the  3000-4000&  solar  uv  radiation. 


SUMMARY 


The  atmospheric  heating  rates  due  to  and  0^  absorption  of 
solar  radiation  have  been  successfully  parameterized  with  an  accuracy 
of  + 5 percent  from  15  to  120  km.  In  addition,  the  diffuse  solar 
radiation  produced  by  multiple  scattering  and  ground  reflection  has 
been  adequately  described  with  a simple  radiative  transfer  model  of  a 
purely  absorbing  layer  on  top  of  a reflecting  layer.  These  parameteri- 
zations  are  suitable  for  use  in  complex  numerical  models  of  atmospheric 
circulation. 
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xd) 

-2  -1 
Flux  (ergs  cm  see  ) 

O(02) (cm2) 

a(03) (cm2 

7500 

1.27(4) 

3.5(-22) 

7400 

1.30(4) 

3.8(-22) 

7300 

1.34(4) 

4. 2 (-22) 

7200 

1.37(4) 

5. 4 (-22) 

7100 

1.40(4) 

6.5(-22) 

7000 

1.43(4) 

8.2(-22) 

6900 

1.47(4) 

1.0(-21) 

6800 

1.50(4) 

1. 3(-21) 

6700 

1.54(4) 

1 . 6(-21) 

6600 

1.54(4) 

2.0(-21) 

6500 

1.56(4) 

2.4(-21) 

6400 

1.62(4) 

2.9 (—2 1 ) 

6300 

1.65(4) 

3.4(-21) 

6200 

1.68(4) 

3.9 (—21 ) 

6100 

1.72(4) 

4. 5 (-21) 

6000 

1.75(4) 

4.6 (—2 1 ) 

5900 

1.78(4) 

4.0 (—2 1 ) 

5800 

1.82(4) 

4.3(-21) 

5700 

1.83(4) 

4.3 (—2 1 ) 

5600 

1.83(4) 

3.5(-21) 

5500 

1.85(4) 

3.1 (—2 1 ) 

5400 

1.89(4) 

2 . 7 (-21) 

5300 

1.89(4) 

2.3(-21) 

5200 

1.85(4) 

1.7 (-21) 

5100 

1.87(4) 

1 • 5 (—2 1 ) 

5000 

1.93(4) 

9 . 4 ( — 22 ) 

4900 

1.92(4) 

6.7 (-22) 

4800 

1.94(4) 

5.7 (-22) 

4700 

1.99(4) 

2.7 (—22) 

4600 

2.00(4) 

2.4(-22) 

4500 

1.97(4) 

1 . 2 (-22) 

3550 

9.5(3) 

3.5 (— 22 ) 

3450 

8.9(3) 

1 . 0 (—2 1 ) 

3350 

9.0(3) 

4. 2 (-21) 

3250 

7.8(3) 

1.5 (-20) 

3150 

6.4(3) 

6.2(-20) 

3075 

2.6(3) 

1.5 (—19) 

3000 

3.9(3) 

3.5 (—19) 

2925 

5.2(3) 

1.2 (—18) 

2825 

2.4(3) 

3.3(-18) 

2725 

2.4(3) 

6.8 (-18) 

2625 

1.9(3) 

1.0(— 17) 

2525 

8.1(2) 

1 .1 (-17) 

2450 

3.5(2) 

9.8(-18) 

2400 

2.9(2) 

5 . 0(— 25) 

8.2 (—18) 

2350 

2.9(2) 

1.4(-24) 

6.4 (—18) 

2300 

3.1(2) 

3.3(-24) 

4.5(-18) 

2250 

3.4(2) 

5.3(-24) 

3.0(-18) 
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Table 

L:  Solar  Fluxes  and  Absorption  Cross  Sections*  - Cont. 

A (ft) 

-2  -1 
Flux  (ergs  cm  sec  ) 

a(02)(cm2) 

a(03) (cm2) 

2200 

2.9(2) 

7.6(-24) 

1.8 (—18) 

2150 

2.4(2) 

9.6(-24) 

1.0 (—18) 

2100 

1.5(2) 

1.1 (-23) 

5.5(-19) 

2060 

2.0(1) 

1.3 (—23) 

3.8 (—19) 

2025 

1.7(1) 

See  Table  2 

3.3(-19) 

2000 

1.4(1) 

tl 

3.1 (—19) 

1985 

7.0(0) 

It 

3.3C-19) 

1972 

5.8(0) 

It 

3 . 6(— 19) 

1947 

1.0(1) 

II 

4.1(-19) 

1924 

7.0(0) 

II 

4.5(-19) 

1902 

6.5(0) 

II 

5 . 2 (—19) 

1882 

5.2(0) 

II 

5 . 9 (—19) 

1863 

4.2(0) 

II 

6.5(-19) 

1846 

2.8(0) 

•1 

6. 8(-19) 

1830 

2.7(0) 

II 

7 . 0(-19) 

1816 

2.5(0) 

tl 

7.3(-19) 

1804 

1.8(0) 

It 

7.5 (—19) 

1793 

1.3(0) 

tl 

7 • 7 (-19) 

1783 

1.2(0) 

tl 

7 . 9 (—19) 

1775 

8.1(-1) 

tl 

8. 1 (-19) 

1769 

6. 7 (-1) 

tl 

8.2 (—19) 

1763 

1.3(0) 

It 

8.2 (—19) 

1750 

1.0(0) 

It 

8.3(-19) 

1740 

1.7(0) 

3.7 (-19) 

1720 

1.5(0) 

5 . 9 (— 1 9) 

1700 

1.3(0) 

8.5(-19) 

1680 

8.7 (-1) 

1.2 (—18) 

1660 

7 . 7 (-1) 

1.8 (—18) 

1640 

5.6 (— 1 ) 

2.5 (—18) 

1620 

4.0(-l) 

3.4(-18) 

1600 

3.2(-l) 

4 . 7 (-18) 

1580 

2.9 (— 1 ) 

6. 0(— 18) 

1560 

3.4(-l) 

7.3  C~ 18 ) 

1540 

3.3 (—1) 

8. 5(-18) 

1520 

2.0(-l) 

1. 0(-17) 

1500 

1.5 (— 1) 

1.1 (-17) 

1480 

1.3 (— 1 ) 

1.2 (—17) 

1460 

l.l(-l) 

1 • 3(— 17) 

1440 

8.3(-2) 

1.5 (—17) 

1420 

7 .0(-2) 

1.5 (—17) 

1400 

1-3 (—1 ) 

1.4 (—17) 

1380 

4. 8 (-2) 

1-3 (—17) 

1360 

6. 3(-2) 

8.0(-18) 

1340 

1.2 (— 1) 

2.3(-18) 

1320 

4.3(-2) 

1.4(-18) 

1300 

7.6 (—2 ) 

5 . 0(-19) 

1280 

6.9(-3) 

2.8(-19) 

1260 

1.6(-2) 

4.3(-19) 

*The  solar  fluxes  are  integrated  over  wavelength  intervals  whose  center 
value  is  given.  Note  that  the  values  in  parentheses  are  the  powers  of 
10  by  which  the  primary  tabular  entry  is  to  be  multipled. 
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These  coefficients  were  supplied  by  McConnell  (1974). 


2 (km) 

T(°k) 

Table  3: 
[0] (cm-3) 

Model  Atmosphere 
[02](cm“3) 

[N2](cm“3 

[03](cm"3) 

15 

216.6 

3(6) 

8.586(17) 

3.191(18) 

2(12) 

20 

216.6 

1(7) 

3.920(17) 

1.457(18) 

3.1(12) 

25 

221.7 

5(7) 

1.720(17) 

6.391(17) 

4.5(12) 

30 

230.7 

2(8) 

7.823(16) 

2.908(17) 

3.1(12) 

35 

241.5 

3.5(8) 

3.651(16) 

1.357(17) 

1.6(12) 

40 

255.3 

1.5(9) 

1.752(16) 

6.513(16) 

6.7(11) 

45 

267.7 

3(9) 

8.794(15) 

3.269(16) 

2.1(11) 

50 

271.6 

9(9) 

4.660(15) 

1.732(16) 

7(10) 

55 

263.9 

2(10) 

2.565(15) 

9.535(15) 

2.2(10) 

60 

249.3 

2.5(10) 

1.414(15) 

5.255(15) 

7(9) 

65 

232.7 

2.3(10) 

7.564(14) 

2.812(15) 

2.2(9) 

70 

216.2 

2(10) 

3.863(14) 

1.43^15) 

7(8) 

75 

205.0 

4(10) 

1.824(14) 

6.792(14) 

2(8) 

80 

195.0 

6.2(10) 

8.204(13) 

3.104(14) 

3.1(8) 

85 

185.1 

1.4(11) 

3.548(13) 

1.365(14) 

1.2(8) 

90 

183.8 

1.66(11) 

1.422(13) 

5.585(13) 

2.7(7) 

95 

190.3 

1.91(11) 

5.483(12) 

2.218(13) 

4.6(6) 

100 

203.5 

4.150(11) 

1.991(12) 

8.710(12) 

1.3(6) 

105 

228.0 

4.436(11) 

6.653(11) 

3.597(12) 

1.5(5) 

110 

265.5 

3.228(11) 

2.500(11) 

] .585(12) 

1.4(4) 

115 

317.1 

2.148(11) 

1.086(11) 

7.447(11) 

1.2(3) 

120 

380.6 

1.422(11) 

5.420(10) 

3.793(11) 

1.5(2) 

Table  4:  Conversion  Factors  for  °K  day  to  ergs  cui  sec" 


Z(km) 

15 

2.26(-2) 

20 

1.03 (-2) 

25 

4.52(-3) 

30 

2.06(-3) 

35 

9. 60 (-4) 

40 

4. 60(-4) 

45 

2.31(-4) 

50 

1 . 23(-4) 

55 

6.75 (—5) 

60 

3. 7 3 (—5) 

65 

1. 99 (—5) 

Z (km) 

70 

1 ,02(-5) 

75 

4.80(-6) 

80 

2. 19 (-6) 

85 

9.59(-7) 

90 

3.90(-7) 

95 

1 - 54 (—7 ) 

100 

6.01 (-8) 

105 

2.43(-8) 

110 

1 . 06(-8) 

115 

5.07 (-9) 

120 

2.63 (—9) 

23 


SUMMER  LATITUDE  (DEGREES)  WINTER 

- Contours  of  constant  estimated  heating  rate  for  solstitial  conditions. 
Model  atmosphere  s assumed  invariant  with  latitude. 
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COMMANDER 

U.  S.  Army  Missile  Command 
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COMMANDER 

U.  S.  Army  Nuclear  Agency 
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ATTN:  John  E . Doncarlos 
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ATTN:  R.  Eastman 
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Washington,  D.  C.  20390 
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Naval  Research  Laboratory 
Washington,  D.  C.  20375 


ATTN: 

HDQ  COMM  DIR  Bruce  Wald 

ATTN: 
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CODE  7701  Jack  D.  Brown 

ATTN: 
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25  copies 
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COMMANDER 
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Strategic  Systems  Project  Office 
Navy  Department 
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ATTN:  NSP-2l4l" 

COMMATOIR 
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ENT  AF3,  Co.,  30912 
ATTN:  ADDA 

Headquarters 

U.  S.  Army  Elect  Warfare  Lab  (ECOM) 

White  Sands  Missile  Range,  NM  88002 
ATTN:  E.  Butterfield 

AF  Cambridge  Rsch  Labs,  AFSC 
L.  G.  Hanscom  Field 
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ATTN:  OPR  James  C.  Ulwick 

ATTN:  OP  John  S.  Garing 
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AF  Weapons  Laboratory,  AFSC 
Kir t land  AFB,  NM  87117 

ATTN:  John  M.  Kamm  SAS 

ATTN:  SUL 

ATTN:  DYT  CAPT.  David  W.  Goetz 

ATTN:  DYT  MAJ.  Don  Mitchell 

AFTAC 

Patrick  AFB,  FI.  32925 

ATTN:  TF  MAJ.  E.  Hines 

ATTN:  TF/CAPT.  Wiley 

ATTN:  TN 

Air  Force  Avionics  Laboratory,  AFSC 
Wright -Patterson  AFB,  Oh.  45433 

ATTN:  AFAL  AVWE  Wade  T.  Hunt 

Assistant  Chief  of  Staff 
Studies  and  Analysis 
Headquarters,  U.  S.  Air  Force 
Washington,  D.  C.  20350 

ATTN:  LTC.  A.  D.  Dayton 


Headquarters 

Electronics  Systems  Division  (AFSC) 

L.  G.  Han scorn  Field 
Bedford.  Ma.  01750 

I Air:::  XRE  LT.  Michaels 

ATTN:  LTC  J.  Morin  CDEF  XRC 

ATT:;:  YSEV  LTC.  David  C.  Sparks 

COMMANDER 

Foreign  Technology  Division,  AFSC 
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ATTN:  TD-BTA  LIBRARY 

! HQ  USAF/RD 

Washington,  D.  C.  20550 
ATTN:  RDQ 

COMMANDER 

Rome  Air  Development  Center,  AFSC 
Griffith  AFB,  N.  Y.  13440 

ATTN:  EMTLD  Doc  Library 

COMMANDER  IN  CHIEF 
Strategic  Air  Command 
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ATTN:  RDPO  LT.  Alan  B.  Merrill 

Los  Alamos  Scientific  Laboratory 

P.  0.  Box  1663 

Los  Alamos , NM  87544 

ATTN:  DOC  CON  for  R.  F.  Taschek 

ATTN:  DOC  CON  for  Milton  Peek 

Sandia  Laboratories 
P.  0.  Box  5800 
Albuquerque,  NM  87115 

ATTN:  DOC  CON  for  A.  Dean  Thronbrough 

ATTN:  DOC  CON  for  W.  D.  Brown 

ATTN:  DOC  CON  for  D.  A.  Dahlgren,  ORG  1722 

University  of  California 
Lawrence  Livermore  Laboratory 
P.  0.  Box  808 
Livermore,  CA  94550 

ATTN:  Tech  Info  Dept  L-3 

Department  of  Commerce 
National  Oceanic  & Atmospheric  Admin. 
Environmental  Research  Laboratories 
Boulder,  CO  30302 

ATTN:  Joseph  h . Pope 

ATTN:  C.  L.  Rufenach 

Department  of  Commerce 
Office  for  Telecommunications 
Institute  for  Telecom  Science 
Boulder,  CO  80302 

ATTN:  Glenn  Falcon 

ATTN:  G.  Reed 

ATTN:  L.  A.  Berry 

Department  of  Transportation 
Transportation  Rsch . System  Center 
Kendall  Square 
Cambridge,  MA  02142 

ATTN:  TER  G)  Harowles 


NASA 

Goddard  Space  Flight  Center 
Greenbelt,  Md  20771 

ATTN:  CODE  750  T.  Golden 


NASA 

600  Independence  Ave.,  S.  W. 
Washington,  D.  C.  20346 
ATTN:  M.  Dub in 

ATTN:  J.  Holtz 

Aerodyne  Research,  Inc. 
Tech/Ops  Building 
20  South  Avenue 
Burlington,  MA  01803 
ATTN : M . Camac 

ATTN:  F.  Bien 

Aerospace  Corporation 
P.  0.  Box  92937 
Los  Angeles,  CA  90009 


ATTN: 

T.  M.  Salmi 

ATTN: 

S . P . Bower 

ATTN: 

V.  Josephson 

ATTN: 

SMFA  for  PWW 

ATTN: 

R.  Grove 

ATTN: 

R.  D.  Rawcliffe 

ATTN: 

T.  Taylor 

ATTN: 

Harris  Mayer 

ATTN: 

D.  C.  Cartwright 

Analytical 

Systems  Corporation 

25  Ray  -Avenue 

Burlington. 

MA  0l90~ 

A" — 'r: 

Radio  Sciences 

Avco-Iveret 

t Research  Laboratory 

23£3  Revere 

Beach  Parkway 

Everett,  MA  0211-9 

ATTN : 

Richard  M.  Patrick 

Boeing  Company,  The 

P.  0.  Box  3707 

Seattle,  WA  93124 

ATTN: 

D . Murray 

ATTN: 

Glen  Keister 

Brown  Engineering  Company,  Inc. 
Cummings  Research  Park 
Huntsville,  AL  35807 

ATTN:  David  Lambert  MS  18 

California  at  San  Diego,  Univ.  of 
Building  500  Mather  Campus 
3172  Miramar  Road 
La  Jolla,  CA  92037 

ATTN:  Henry  G.  Booker 
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Calspan 
P.  0.  Box  235 
Buffalo,  N.  Y.  14221 

ATTN:  Romeo  A.  Deliberis 

Computer  Sciences  Corporation 

P.  0.  Box  530 

6565  Arlington  Blvd. 

Falls  Church,  VA  22046 
ATTN:  H.  Blank 

ATTN:  Barbara  F.  Adams 

Comsat  Laboratories 
P.  0.  Box  115 
Clarksburg,  Md . 20734 

ATTN:  R.  R.  Taur 

Cornell  University 

Department  of  Electrical  Engineering 

Ithaca,  N.  Y.  14850 

ATTN:  D.  T.  Farley,  Jr. 

ESL,  Inc. 

495  Java  Drive 
Sunnyvale,  CA  93102 
ATTN:  J.  Roberts 

ATTN:  V.  L.  Mower 

ATTN:  James  Marshall 

ATTN:  R.  K.  Stevens 

General  Electric  Corn any 
Tempo-Center  for  Advanced  Studies 
Sic  State  Street 
Santa  3arbara,  CA  35102 
ATTN:  Don  Chandler 

ATTN:  DASIAC 

ATTN:  DASIAC  Art  Feryok 

ATTN : Warren  S . Knapp 

General  Electric  Company 
P.  0.  Box  1122 
Syracuse,  N.  Y.  13201 

ATTN:  F.  A.  Reibert 

General  Research  Corporation 
P.  0.  Box  3587 
Santa  Barbara,  CA  93405 
ATTN:  John  Ise,  Jr. 


Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 

ATTN : Technical  Library 

ATTN:  Neil  Brown 

ATTN:  T.  N.  Davis 

GTE  Sylvania,  Inc. 

I89  B Street 

Needham  Heights,  MA  02194 
ATTN:  Marshall  Cross 

HRB -SINGER,  Inc. 

Science  Park,  Science  Park  Road 

P.  0.  Box  60 

State  College,  PA  16801 

ATTN : Larry  Feathers 

Honeywell  Incorporated 
Radiation  Center 
2 Forbes  Road 
Lexington,  MA  02173 

ATTN:  W.  Williamson 

Illinois,  University  of 
Department  of  Electrical  Engineering 
Urbane,  XL  6 1001 

AITS:  K.  C.  Yah 

Institute  for  Defense  Analyses 
1*00  Army-Navy  Drive 
Arlington,  VA  22202 

ATTIC:  Ernest  Bauer 

ATTN:  Hans  Wolfhard 

ATTN:  J.  M.  Aein 

ATTIC:  Joel  Bengston 

Inti  Tel  St  Telegraph  Corporation 
500  Washington  Avenue 
Nutley,  N.  J.  07 HO 

ATTN:  Technical  Library 

ITT  Electro-Physics  Laboratories,  Inc. 
911*0  Old  Annapolis  Road 
Columbus,  Md.  21043 

ATTN:  John  M.  Kelso 
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Johns  Hopkins  University 
Applied  Physics  Laboratory 
8621  Georgia  Avenue 
Silver  Spring,  MD  20910 

ATTN:  Document  Lib-^^ian 

Lockheed  Missiles  & Space  Co.,  Inc. 

P.  0.  Box  504 
Sunnyvale,  CA  94388 
ATTN:  Dept.  60- 12 

Lockheed  Missiles  and  Space  Company 
3251  Hanover  Street 
Palo  Alto,  CA  94304 

ATTN:  Billy  M.  McCormac,  Dept  52-14 

ATTN:  Martin  Walt,  Dept  52-10 

ATTN:  Richard  G.  Johnson,  Dept  52-12 

ATTN:  J.  B.  Reagan,  Dept  52-12 

ATTN:  John  Kumer 

ATTN:  Robert  D.  Sears,  Dept  52-14 

MIT  Lincoln  Laboratory 
P.  0.  3ox  73 
Lexington,  MA  02173 

ATTN:  Mr.  Walden,  X113 

ATTN:  D.  Clark 

ATTN:  James  H.  Pannell,  L-2 U6 

ATTN : Lib  A-032  for  David  M.  Towle 

Martin  Marietta  Corporation 
Denver  Distribution 
P.  0.  Box  179 
Denver,  CO  30201 

ATTN:  Special  Projects  Program  248 

Maxwell  Laboratories,  Inc. 

9244  Balboa  Avenue 
San  Diego,  CA  92123 

ATTN:  A.  J.  Shannon 

ATTN : V . Fargo 

ATTN:  A.  N.  Rostocker 

McDonnell  Douglas  Corporation 
5301  Bo Isa  Avenue 
Huntington  Beach,  CA  92657 
ATTN : J . Moule 

ATTN:  N.  Harris 
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Physical  Dynamics,  Inc. 

P.  0.  Box  1069 
Berkeley,  CA  94701 

ATTN:  Joseph  B.  Workman 

Physical  Sciences,  Inc. 

607  North  Avenue,  Door  18 
Wakefield,  MA  01880 
ATTN : Kurt  Wray 

R & D Associates 

P.  0.  Box  358O 

Santa  Monica,  CA  90433 

ATTN:  Robert  E.  Lelevier 

ATTN:  Forest  Gilmore 

ATTN:  Richard  Latter 

ATTN:  William  B.  Wright,  Jr. 

R & D Associates 
1815  N.  Ft.  Myer  Drive 
11th  Floor 
Arlington,  VA  2209 

ATTN:  Herbert  J.  Mitchell 

Rand  Corporation,  The 
1700  Main  Street 
Santa  Monica,  CA  90a05 

ATTN:  Cullen  Crain 

ATT.' : James  Oakley 


Science  Applications.  Inc. 

?.  C.  3cx  23pl 
La  Jolla,  CA  92035 

ATT.' : Daniel  A.  Hamlin 

ATTN:  D.  Sachs 

ATTN:  E.  A.  Straker 

Space  Data  Corporation 
1331  South  26th  Street 
Phoenix,  AZ  35034 

ATTN:  Edward  F.  Allen 


Stanford  Research  Institute 
333  Ravenswood  Avenue 
Menlo  Park,  CA  9^*025 
ATTN:  M.  Baron 

ATTN:  L.  L.  Cobb 

ATTN:  Walter  G.  Chestnut 

ATTN:  David  A.  Johnson 

ATTN:  Charles  L.  Rino 

ATTN:  E.  J.  Fremouw 

ATTN:  Ray  L.  Leadabrand 

Stanford  Research  Institute 
306  Wyyn  Drive,  N.  W. 

Huntsville,  AL  35805 

ATTN:  Dale  H.  Davis 

Technology  International  Corporation 
75  Wiggins  Avenue 
Bedford,  MA  01730 

ATTN:  W.  P.  Boquist 

TRW  Systems  Group 
One  Space  Park 
Redondo  Beach,  CA  91273 
.-.777 : P.  H.  Katsos 

-777 : j.  w.  Levery 


Utah  Stpts 

Universitv 

Logan,  UT 

31321 

att::: 

C.  Wyatt 

ait::: 

D.  Burt 

ATTN: 

Kay  Baker 

ATTN: 

Doran  Baker 

Visidyne, 

Inc. 

19  Third  Avenue 

North  West 

Industrial  Park 

Burlington 

, MA  01803 

ATTN: 

William  Reidy 

ATTN: 

Oscar  Manley 

ATTN: 

T.  C.  Degges 

ATTN: 

J.  W.  Carpenter 
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